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Abstract. We analysed the XMM-Newton archival observations of 16 neutron star (NS) low-mass X-ray binaries 
(LMXBs) to study the Fe K emission in these objects. The sample includes all the observations of NS LMXBs 
performed in EPIC pn Timing mode with XMM-Newton publicly available until September 30, 2009. We performed 
a detailed data analysis considering pile-up and background effects. The properties of the iron lines differed from 
previous published analyses due to either incorrect pile-up corrections or different continuum parameterization. 
80% of the observations for which a spectrum can be extracted showed significant Fe line emission. We found 
an average line centroid of 6.67 ±0.02 keV and a finite width, a, of 0.33 ±0.02 keV. The equivalent width of 
the lines varied between 17 and 189 eV, with an average weighted value of 42 ± 3 eV. For sources where several 
observations were available the Fe line parameters changed between observations whenever the continuum changed 
significantly. The line parameters did not show any correlation with luminosity. Most important, we could fit the 
Fe line with a simple Gaussian component for all the sources. The lines did not show the asymmetric profiles that 
were interpreted as an indication of relativistic effects in previous analyses of these LMXBs. 

Key words. X-rays: binaries - Accretion, accretion disks - X-rays: individual: 4U 1608— 52, CenX— 4, 
AqlX-1, IGRJ00291+5934, 4U 1705-44, XTE J1807-294, 4U 1728-34, 4U 1543-62, 4U 0614+09, 4U 1636-536, 
SAXJ1808.4-3658, GX340±0, SerX-1, 4U 1735-44, GX9±9, GX 349+2 



1. Introduction 

Accreting binaries often show iron line emission in their 
X-ray spectra. The ability of Chandra, XMM-Newton, 
Suzaku and Swift to obtain medium to high resolution 
spectra covering the iron K energy band has opened a new 
era in the observations of stellar-mass black holes (BHs) 
and NS binaries. The large effective area of these observa- 
tories is crucial for the detection of iron line emission while 
the high resolution enables us to discriminate between 
narrow and intrinsically broad features. Br oad iron lines 



sponsible for it are still controversial. Undoubtedly the 
most exciting possibility is that such lines originate in 
the disc close to the BH event horizon or to the NS by 
fluorescent emission following illumination (and photoin- 
isatio n) of the accretion disc by an external source of X- 
ravs (iRevnolds fc Nowa^ bOOa iFabian fc Miniuttil 12005 



Mattll2006r ). In this model, the combination of relativistic 



are a l so present in Active Galactic Nuclei ([Tanaka et al. 



1995t iNandra et al.l Il997l ) though studies of large sam- 



ples of sources show that the ave rage fraction of broad 
Fe lines is never hi gher than 40% ( Guainazzi et al.ll2006t 



Nandra et al.ll2007 ) 



The origin of broad iron li nes has been extensively dis- 



cussed in the literature (e.g. ISunvaev fc Titarchuk 



Fabian et al."l989^' Done et al.ll2007l: IRoss fc FabianI 



198C; 



2007 



Titarchu k ct al. 200!|). However, even after the advent of 



the current powerful X-ray observatories, the exact de- 
termination of the line width and the mechanisms re- 



Doppler effects arising from the high orbital velocities and 
gravitational effects due to the strong gravitational field 
in the vicinity of the compact object smear the reflected 
features. In this case, detailed X-ray spectroscopy of iron 
line features can be used to study Doppler and gravita- 
tional redshifts, thereby providing key information on the 
location and kinematics of the material in the vicinity of 
the compact object. Most interesting of all is the potential 
for establishing BH spin using relativistic iron Ka lines. 
The spin value is constrained mainly by the lower bound- 
ary of the broad line, which depends on the inner bound- 
ary of the disc emission (identified with the marginally 
stable orbit) where the gravitational redshift is maximal. 
Therefore, lines which are strongly skewed toward lower 
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energies can indicate black hole spin ( Miller et al.ll2009t 
Hiemstra et"aLll2009HReis et al.ll20Q9t ). 

The recent claim of broad skewed iron lines from NS 
binaries has opened the exciting possibility of determining 
an upper limit to the radius of the NS, the most difficult 
parameter to o btain in order to constrain the equ ation 
of state of NSs (JBhattacharvva fc Strohmaveiil2007l) . The 
advantage of using iron K emission lines as a probe of the 
NS radius is that they only require short observations to 
clearly reveal the relativistic lines, and do not require any 
knowledge of the distance to the object. 

An alternative location of the line emitting region 
could be the inner part of the so-called Accretion Disc 
Corona (ADC), formed by evaporation of the outer lay- 
ers of the di sc illuminated by the emission of the central 
object (e.g., IWhite fc Holtl Il982l) . In the corona, where 
the gas is highly ionised, the iron line emission is likely 
to come from recombination onto hydrogenic and helium- 
like iron. In the region below the transition region the 
iron will be at most a few ti mes ionised and fluoresce nce 
will be the dominant process ( Kallman fc Whitelll989r) . In 
the corona, Compton scattering will broaden and shift the 
centroid of the line. The line broadening and shift are due 
to Doppler and/or recoil effects, depending on the tem- 
perature of the plasma c ompared to the energy of the 
photo ns (jPozdnvakov et al. 1979; ISunvaev &: TitarchukI 
19801) . The final line width is determined by the simul- 
taneous effects of blending, Compton scattering and ro- 
tation. For a value of /~0.1 (/ being the fraction of the 
incident X-ray flux which is ass umed to penetrate t h e base 
of the hot medium or corona), iKallman fc White! (|l989r) 
calculated a centroid energy of 6.6-6.7 keV for the iron 
line and a width of ^0.6 keV at small radii (R ^ 10® cm). 
For /~ 1, the width can be as large as ~1 keV (FWHM). 
Alternatively, widths of ~1 keV can also be attained via 
emission from a photoionised gas with large Thomson 
depth (r ,^3). Such large Compton depths could arise 
if the corona is heated by some mechanism other than 
X-rays, or if the line of sight of the observer through the 
corona were very diffe rent from that of the compact X-ray 
source (JKallman fc W hite 1989). 

A third scenario for the formation of broad iron 
lines in LMXBs could be one where extensive red wings 
form by recoil of line photons in an optically thick 
medi um expanding, or converging, at relativistic veloc- 
ities dri tarchukct al l 120031 iLaming fc TitarchukI 12004 : 
Laurent fc Titarchuk l2007t ). It is assumed that the wide- 
open wind is launched at some disc radius where presum- 
ably the local radiation force exceeds the local disc gravity. 
The wind should be illuminated by the emission of X-rays 
formed in the innermost part of the source. The Ka line 
is generated in a narrow wind shell and the line profile 
is formed in the partly ionised wind when Ka line pho- 
tons are scattered off the diverging flow (wind) electrons 
( Titarchuk et al.ll2009t ). The red-skewed part of the spec- 
trum is formed by photons undergoing multiple scatterings 
while the primary peak is formed by photons escaping di- 
rectly to the observer. 



The line properties of an ADC are in general ex- 
pected to be differeii t from those of an accretion disc 
( Brandt fc Mattlll994f ). In the former case the equivalent 
width (EW) is nearl y independent of in clination except 



at very high angles (jVrtilek et al.l Il993r ) , while the de 



pendence on inclination of the case of an accretion disc 
is significant (JBrandt fc Mattlll994l) . Therefore iron lines 
appear to be a powerful tool for distinguishing between 
an ADC, an accretion disc, or any other origin of the iron 
line in LMXBs. In any case, a proper modelling of the iron 
line is mandatory to determine uniquely the origin of the 
line and, in the case of relativistic broadened iron lines, to 
derive the spin of the BH or an upper limit to the radius 

of the NS. 

White et al.l ( 19861) performed a systematic analysis of 



six NS LMXBs with EXOSAT and found line emission in 
83% of the objects. The hues were broad with FWHM 
of 0.8-1.3 keV and had centroid energies of 6.6-6.9 keV. 
Further the line properties di d not show o bvious correla- 
tion with luminosity. This led IWhite et al.l (|l986: ) to con- 
clude that the only plausible broadening mechanism was 
Comptonization in a cloud with a Thomson depth of a few 



and a n electron temperature close to 7 kcV. Hir ano et al. 



(|19871 ) analysed ten NS LMXBs observed with Tenma and 
found line emission from six of them. The line energy had a 
weighted average of 6.66 ±0.05 keV with EWs in the range 
20-60 eV. A significant line width of 0.55 keV was only ob- 
tained for one source, though a width of '^l keV could not 
be rejected for the other cases. They attributed the line 
emission to recombination processes in the accretion disc 
corona and explained its width as Compoton scattering in 
the region with a n, optical depth of 1 and an average tem- 
perature ~1 keV. lAsai et al.l (jgOOOJ) extended the analysis 
of the iron K emission lines to twenty NS LMXBs ob- 
served with ASCA. They detected significant iron lines 
from about half of the sources. The average properties of 
the lines were a line centre of 6.56 keV and a finite width of 
~0.5 keV (FWHM) in six of the sources. They also found 
a large scatter in the EW of the detected lines, ranging 
between 10 and 170 eV. They concluded that the iron K 
lines are likely produced through the radiative recombina- 
tion of photoionised plasma and explained the line width 
as a combination of line blending, Doppler broadening and 
Compton scatterings. 

The picture outlined by the analyses based on 
EXOSAT, Tenma and ASCA data changed with re- 
cent analyses of observations of NS LMXB obser- 
vations by XMM-Newton, Chandra and Suzaku. A 
number of skewed iron lines have been claimed and 
their origin attributed to relat i vistic effects (e.g. . 



Bhattacharvva fc Strohmav er 2007; IPaoitto et al.l 120091 : 



Di Salvo et al.l2009HD'Aret al.ii20090 or to Compton scat 



tering in an optically thick medium e xpanding, or converg- 



ing, at relativistic velocitie s (e.g., iLaurent fc Titarchuk 



20071: iTitarchuket al.ll2009l ). However only some of the 



NS LMXBs showed such as ymmetric lines whil e others 
showed symmetric lines (^e.g.. lCackett et al.ll2008l ). A sys- 
tematic analysis of the iron K lines is at this stage fun- 
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damental in order to establish why some sources exhibit 
skewed lines while others show symmetric lines or why 
lines are only present in some of the sources. Having an 
answer to these questions would help to clarify the cur- 
rently controversial origin of the lines. 

With this aim we performed a systematic analysis of 
all the NS LMXBs observed by XMM-Newton since the 
beginning of the mission and publicly available up to the 
30th of September 2009. We excluded from this sample 
observations of dipping or ADC sources, i.e. with known 
inclinations above 70°, since their analyses at the Fe band 
are complicated due to strong absorption in the line of 
sight. We chose XMM-Newton to perform this study for 
the following reasons: firstly, its high effective area both 
below and above the iron K band allows a good determi- 
nation of the shape and width of the iron line as well as the 
continuum. Secondly, XMM-Newton has provided most of 
the detections and well determined profiles of broad iron 
lines. Finally, XMM-Newton, after 10 years of operation, 
has an ample archive of public data which can be used to 
determine class properties in the least unbiased way. 

In total we analysed 26 observations of 16 sources from 
the XMM-Newton archive. We restricted our analysis to 
sources observed in the EPIC pn Timing Mode. This mode 
is specially suited for the observation of bright sources. 
Therefore, by selecting all the observations performed in 
this mode we obtained a sample of spectra with the best 
possible statistics. 

In what follows we first present the properties of the se- 
lected sample and our analysis method. Special attention 
was payed to eliminate pile-up effects and to the treat- 
ment of background. Then, we performed spectral fitting 
accounting for the excess emission at the Fe K energy 
band. We finally discuss the implications of our analysis 
regarding the characteristics of the lines and their possible 
origin. In an appendix we compare the characteristics of 
the lines in this work with previous analyses of the same 
observations and discuss the reason for the discrepancies 
whenever found. 

2. Observations and data reduction 



The XMM-Newton Observatory (jJansen et al.l 120011 ) in 



eludes three 1500 cm^ X-ray telescopes each with an 
EPIC (0.1-15 keV) at the focus. Tw o of the EPIC ima^ 
ing spectrometers use MPS CCDs d 
and one uses p n CCDs dStriider et al. 



(0.35-2.5 keV, iDen Herder et al.l 1200 



urner et al.l 120011) 



20011) 



J) 



The RGSs 
are located be- 
hind two of the telescopes. In addition, there is a co- 
aligne d 30 cm diameter Optical/UV Monitor telescope 
(OM, iMason et al.l 120011 ). providing simultaneously cov- 
erage with the X-ray instruments. Data products were re- 
duced using the Science Analysis Software (SAS) version 
9.0. The EPIC MOS cameras were not used in most of the 
observations analysed in this paper because their teleme- 
try was allocated to the EPIC pn camera to avoid Full 
Scientific Buffer in the latter. Since the pn has an effec- 
tive area ~5 times higher at 7 keV than the MOS CCDs 



and the latter were not available for most of the obser- 
vations, we present in this work only the analysis of the 
EPIC pn data. We did not analyse the RGS data since 
they do not cover the Fe K energy band in which we are 
interested. 

Table [T] is a summary of the XMM-Newton observa- 
tions. The EPIC pn was used in Timing Mode for all the 
observations. In this mode only one CCD chip is oper- 
ated and the data are collapsed into a one-dimensional 
row (4f4) and read out at high speed, the second dimen- 
sion being replaced by timing information. This allows a 
time resolution of 30 fis. We used the SAS task epf ast on 
the event files to correct for a Charge Transfer Inefficiency 
(CTI) effect which has been seen in the EPIC pn Timing 
mode when high count rates are present Ij. Ancillary re- 
sponse files were generated using the SAS task arfgen 
following the recommendations of the XMM-Newton SAS 
User guide |3 for piled-up observations in Timing mode 
whenever applicable. Response matrices were generated 
using the SAS task rmf gen. We extracted one EPIC pn 
spectrum per observation, not taking into account any 
intra-observational variability (see Sect. 13.1) ). Bursts were 
excluded for the calculation of the total energy spectra 
whenever present. 

Light curves were generated with the SAS task 
epiclccorr, which accounts for time dependent correc- 
tions within a exposure, like dead time and GTIs. 

2.1. Pile-up treatment 

Since the average count rate in the EPIC pn was close to, 
or above, the 800 counts s~^ level at which pile-up effects 
may become significant for at least eight observations in 
the sample, we investigated in detail the presence of pile- 
up before extracting the spectra. We used the SAS task 
epatplot, which utilizes the relative ratios of single- and 
double-pixel events deviating from the standard values as 
a diagnostic tool in case of significant pile-up in the pn 
camera Timing mode. We found that the spectra of twelve 
observations from eight sources were affected by pile-up 
(see Table [1]). 

Then, we extracted several spectra selecting single and 
double Timing Mode events (patterns to 4) and dif- 
ferent spatial regions for each of the piled-up observa- 
tions. For 4U 1705-44 (Obs 0551270201), 4U 1636-536, 
4U 1735-44 and GX 9-h9 (Obs 0090340601) the source co- 
ordinates fell in the centre of the central column of the 
CCD. Therefore, source events were first extracted from a 
64" (15 columns) wide box centred on the source position 
(Region 1). Then, we excluded the neighbouring 1, 3, 5 
and 7 columns from the centre of Region 1 (Regions 2- 



^ More information about the CTI correction can be found 
in the EPIC status of calibration and data analysis and in the 
Current Calibration File (CCF) release note Rate- dependent 
CTI correction for EPIC-pn Timing Modes, by Guainazzi et al. 
(2008), at http://xmm.esac.esa.int/external/xmm_calibration 

^ http://xmm.esac.esa.int 
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Table 1. XMM-Newton archival observations of NS LMXBs performed in EPIC pn Timing mode ordered by Right 
Ascension. T is the total effective EPIC pn exposure time. C is the pn 0.7-10 keV persistent emission count rate 
after dead time correction, calculated as the mean of the Gaussian function used to fit the count rate distribution (cr 
represents the standard deviation of the distribution). The source class is A (atoll), Z (Z source), UCXB (ultra-compact 
X-ray binary) or AMXP (accreting millisecond X-ray pulsar). T means transient system. The last column shows the 
number of columns excised from the centre of the PSF to remove the pile-up effects in spectral analysis. 





Source 


Class 


Observation 


Observation Times (UTC) 


T 


C[o] 


Pile- up 


Columns 






ID 


Start 
(year mon day) 


(hr:mn) 


End 
(hr:mn) 


(ks) 


(s-^) 




removed 


4U 0614-^09 


A,UCXB 


0111040101 


2001 Mar 13 


12:27 


17:11 


10 


240 [8] 


N 





CenX-4 


T 


0144900101 


2003 Mar 01 


15:11 


14:48 


- 


- 


- 


- 


4U 1543-62 


UCXB 


0061140201 


2001 Feb 04 


13:16 


03:13 


46 


203 [5] 


N 





4U 1608-52 


A, T 


0074140101 


2002 Feb 13 


16:04 


20:57 


- 


- 


- 


- 






0074140201 


2002 Feb 15 


01:34 


06:22 


- 


- 


- 


- 


4U 1636-536 


A 


0303250201 


2005 Aug 29 


17:47 


02:47 


29 


243 [15] 


N 









0500350301 


2007 Sep 28 


15:07 


00:17 


19 


507 [11] 


Y 


1 






0500350401 


2008 Feb 27 


03:38 


15:01 


37 


652 [24] 


Y 


3 


GX 340+0 


Z 


0505950101 


2007 Sep 02 


12:40 


02:34 


40 


801 [50] 


Y 


8 


GX349-f2 


Z 


0506110101 


2008 Mar 19 


15:07 


23:00 


7 


2043 [69] 


Y 


8 


4U 1705-44 


A 


0402300201 


2006 Aug 26 


04:27 


14:57 


35 


30 [1] 


N 









0551270201 


2008 Aug 24 


01:41 


17:15 


45 


742 [31] 


Y 


7 


GX9+9 


A 


0090340101 


2001 Sep 04 


09:37 


15:17 


1 


1380 [11] 


Y 


4 






0090340601 


2002 Sep 25 


09:15 


16:10 


5 


1458 [25] 


Y 


5 


4U 1728-34 


A 


0149810101 


2002 Oct 03 


21:48 


05:55 


26 


88 [1] 


N 





4U 1735-44 


A 


0090340201 


2001 Sep 03 


02:57 


09:06 


5 


1198 [23] 


Y 


3 


SerX-1 


A 


0084020401 


2004 Mar 22 


14:58 


21:18 


6 


1074 [23] 


Y 


4 






0084020501 


2004 Mar 24 


14:47 


21:10 


7 


925 [14] 


Y 


4 






0084020601 


2004 Mar 26 


14:18 


20:41 


5 


1014 [32] 


Y 


4 


AqlX-1 


A, T 


0112440101 


2002 Oct 27 


01:03 


03:30 


- 


- 


- 


- 






0112440301 


2002 Oct 15 


01:55 


04:17 


- 


- 


- 


- 






0112440401 


2002 Oct 17 


01:42 


05:52 


- 


- 


- 


- 






0303220201 


2005 Apr 07 


14:30 


18:58 


3 


228 [6] 


N 





IGRJ00291+5934 


AMXP, T 


0560180201 


2008 Aug 25 


04:45 


14:25 


- 


- 


- 


- 


XTEJ1807-294 


AMXP, T 


0157960101 


2003 Mar 22 


13:40 


18:40 


9 


41 [1] 


N 





SAXJ1808.4-3658 


AMXP, T 


0560180601 


2008 Sep 30 


23:15 


17:19 


43 


550 [10] 


Y" 


2 



"For SAX J1808.4— 3658, we found indications for very small pile-up in the two central columns. Therefore we present the 
analysis with both the full PSF and after removal of the 2 central columns in the following sections. 



5, respectively) and extracted one spectrum for each of 
the defined regions. For GX340-hO, GX 349-^2, SerX-1 
and GX 9-^9 (Obs 0090340101) the source feU between two 
columns in the CCD. Therefore, source events were first 
extracted from a 68" (16 columns) wide box centred on 
the source position (Region 1) and Regions 2-5 were de- 
fined by excluding 2, 4, 6, and 8 columns respectively from 
the centre of Region 1. This was done to preserve the best 
symmetry when excluding piled-up events. Table [T] lists 
the number of columns that had to be extracted from each 
source in order to obtain spectra free of pile-up. It is evi- 
dent that pile-up starts to be important already at a count 
rate of >450 s~^. Pile-up depends on the spectral shape 
and the time variation of the source in a complex way. 
Therefore, although the average count rate of an obser- 
vation generally gives an indication of whether the effects 
of pile-up are important, it is also of outmost importance 



to carefully inspect the epatplot to evaluate the PSF ra- 
dius at which the relative ratios of single- and double- pixel 
events do not deviate from the standard values in the full 
energy band. 

As an example of the use of epatplot to determine 
the amount of pile-up, we show the epatplots in a highly 
piled-up source, GX 349-1-2, and in a dim source free of 
pile-up, 4U 1728-34, in Appendix El 

2.2. Background treatment 

In the EPIC pn Timing mode, there is no source-free back- 
ground region, since the PSF of the telescope extends fur- 
ther than the central CCD boundaries. The central CCD 
has a field of view of 13f6 x 4.'4 in the pn. In timing 
mode, the largest column is the one in which the data are 
collapsed into one-dimensional row. Therefore, the maxi- 
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mum angle for backgromid extractfon is 2', com.pared to 
5' for imaging modes. In our sample, sixteen out of the 
nineteen observations (for wliich source emission is de- 
tected) are very bright, with total count rates >200 s^^ 
(see Table [T]). Therefore, the spectra from these sources 
will not be significantly modified by the "real" background 
which contributes less than 1% to the total count rate in 
most of the bandwidth. Conversely, subtracting the back- 
ground extracted from the outer columns of the central 
CCD will modify the source spectrum, since the PSF is 
energy dependent and the source photons scattered to the 
outer columns do not show the same energy dependence 
as the photons focused on the inner columns. 

For the reasons mentioned above, in this work we chose 
not to subtract the "background" extracted from the outer 
regions of the central CCD. This is an appropriate method 
for all the sources for which the "real" background is neg- 
ligible compared to the source count rate at all energies. 
In contrast, for sources with A^h <; 1 x 10^^ atom cm~^, 
the background is expected to contribute significantly at 
energies ^1 keV, where most of the source photons are 
absorbed in the interstellar medium. Therefore for these 
sources it is mandatory to remove the bins where the back- 
ground is significant before spectral analysis. 

With this aim, we extracted the background from the 
outer regions of the central CCD of "blank fields" , where 
we took as examples of "blank fields" the observations of 
AqlX-1 and 4U 1608-52 listed in Table [D for which the 
source was not significantly detected. Then, we inspected 
the spectra of sources with iVn <; 0.5 x 10^^ atom cm^^, 
and compared the "blank field" background with the 
background extracted from the observation for which the 
source is being analysed. We scaled the former so that it 
was comparable to the latter at energies where we did not 
expect a significant flux from the source. We estimated 
such energy as the one at which simulated spectra with 
different values of TVh started to flatten as a consequence 
of being dominated by background events. Then we re- 
moved energy bins where we expect a contribution from 
the background to the total count rate of more than 1%. 
For 4U 1705-44, 4U 1728-34 and GX 3404-0 we removed 
bins below 1.5, 1.8 and 2.2 keV, respectively. We used all 
the other spectra from 0.7 keV up to 10 keV. 

To illustrate the effects of background subtraction in 
the spectral fitting and evaluate the validity of our method 
we chose the weakest and brightest sources of our sample, 
4U 1705-44 (Obs 0402300201) and GX 349-^2. We exam- 
ined the differences in the spectra after (1) not subtracting 
any background from the source spectrum (Spec 1), (2) 
subtracting the background extracted from a blank field 
(Spec 2), and (3) subtracting the background extracted 
from the outer columns of the central CCD during the ob- 
servation (Spec 3). Fig. [T] (left panels) shows the ratio of 
Spec 1 to its best-fit continuum model (black) and the dif- 
ference from Spec 2 (red) and Spec 3 (green) with respect 
to the same model. The red and black residuals are consis- 
tent within the errors at all energies for GX 349-1-2 (lower 
panels) and above ~1 keV for 4U 1705—44 (upper panels). 



The deviation below ~1 keV for 4U 1705—44 is expected, 
since at such energies the background photons have a sig- 
nificant contribution to the total spectrum. In contrast, we 
observe an energy dependent discrepancy of Spec 3 with 
respect to Spec 1 and 2, as expected due to the energy 
dependence of the PSF. These results are consistent with 
the shape of the background shown in Fig. [T] (right pan- 
els) . The "background" extracted from the outer columns 
of the CCD (green) is clearly contaminated by the source 
at all energies (compared to the source spectrum in black), 
while the "background" extracted from a blank field (red) 
shows the shape expected when compared to Fig. 35 of 
the XMM-Newton Users Handbook. 



3. Light curves and spectral fitting 

3.1. Light curves 

Fig. [5] shows EPIC pn light curves and hardness ra- 
tios of all the XMM-Newton observations analysed in 
this paper with a binning of 64 s. The hardness ratio is 
counts in the 3-10 keV energy range divided by those be- 
tween 0.7-3 keV, except for 4U 1705-44, 4U 1728-34 and 
GX340-hO, for which the soft band is 1.5-3 keV, 1.8-3 keV 
and 2.2-3 keV, respectively. The light curves of observa- 
tions with count rates ^ 200 counts s^^ suffer from regu- 
lar telemetry drops which are seen when the light curves 
are plotted with a high time resolution, e.g. 1 s. The SAS 
task epiclccorr accounts for such gaps, as well as for 
dead time with a different origin, by taking into account 
the count rate in adjacent frames (see Sect. [5]). 

The average count rate changed by more than 2 orders 
of magnitude among the studied sources. Observations of 
the same source taken within months showed strong vari- 
ability in some cases (e.g. 4U 1636—536) and a steady level 
in other cases (e.g. ScrX— 1). 

Strong variability of ^30% within one observation 
is only present in the light curves of GX 340-1-0 and 
GX 349-1-2. The variab i hties w ithi n such observation s were 
studied bv lD'Aiet aTl (|2009t ) and llaria etlll (|2009l) . who 
found that the spectral fits to different intervals within the 
observations gave consistent results for the parameters of 
the iron line. 



3.2. X-ray spectra 

We rebinned the EPIC pn spectra to over-sample the 
FWHM of the energy resolution by a factor of 3 and 
to have a minimum of 25 counts per bin, to allow the use 
of the x^ statistic. To account for systematic effects in 
the EPIC pn timing mode we added quadratically a 2% 
uncertainty to each pn spec tral bin. We p erformed spec- 
tral analysis using XSPEC ( Arnaudlll996 ) v ersion 12.3.1. 



We us ed the photo-electric cross-sections of I Wilms et al. 



( 20001 ) to account for absorption by neutral gas with solar 
abundances (the so-called tbabs XSPEC model). Spectral 
uncertainties are given at 90% confidence (Ax^= 2.71), 
and upper limits at 95% confidence. 
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Fig. 1. Left: Ratio of Spec 1 (see Sect. 12. 2p to its best-fit continuum model (black) and the difference from Spec 2 
(red) and Spec 3 (green) with respect to the same model. Spec 3 has been scaled so that the energy bin at 10 keV 
coincides with Spec 1. The original Spec 3 is shown in grey. Right: Source spectrum (black), background from a blank 
field (red) and background extracted from the outer columns of the central CCD (green). The upper and lower panels 
show the results for 4U 1705-44 and GX 349+2 respectively. 



We first fitted the EPIC pn spectra with a model con- 
sisting of a blackbody and a disc blackbody, both modified 
by photo-electric absorption from neutral material (model 
tbabs*(bbodyrad+diskbb) in XSPEC). 

The fits were unacceptable in most of the cases, mainly 
due to the presence of a broad emission feature at ~6- 
7 keV and a significant excess in the pn spectrum be- 
low ^1.5 keV. This excess consisted of a "Gaussian- 
like" component centred at ^1 keV. This feature has 
been previously modelled in a number of sources either 
as an emissio n line or as an e dge, and its nature is 



unclear 



[e:S 



Boirin et al. 



ISidoh et all I2OOII: iBoirin fc Parmaij l2003t 
2OO4L 12005 ). If the feature has an astrophysi- 



cal origin, its energy is consistent with Ne X or a blend of 
Fexx-Fexxiv emission. When the feature is edge- like its 
energy is consistent with O viii. We therefore fitted such 
a feature with a Gaussian component or an edge when- 
ever present and discuss its origin in Sect. [S] We fitted 
the excess at ~6-7 keV with either a Gaussian or a laor 
component (see below). Finally, we added two additional 



Gaussian absorption/emission features at ^^1.84 keV and 
^^2.28 keV. Such features are probably due to an incor- 
rect modelling of the Si and Au absorption in the CCD 
detectors by the EPIC pn calibration and are therefore 
not further discussed. 

Summarising, our final model consisted of disc black- 
body and blackbody components, one Gaussian emission 
feature at ~1 keV (or absorption edge at ~0.87 keV), and 
one emission feature at ~6.5 keV (modelled with Gaussian 
or laor), both modified by photo-electric absorption from 
neutral material, and two narrow Gaussian features at 
^^1.84 and ~2.28 keV to account for the calibration defi- 
ciencies at these energies (Models la and 2a, see Tabled]). 



The fits with Model la were acceptable for fifteen out 
of the nineteen observations for which we extracted a spec- 
trum, with xt between 0.7 and 1.2 for a number of degrees 
of freedom (d.o.f.) between 161 and 221. The parame- 
ters of the best-fit with this model are given in Table [3] 
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Fig. 2. EPIC pn light curves and hardness ratios (HR) for all the observations of our sample with a binning of 64 s. 
Time is shown in hours since the beginning of the observation. The hardness ratio is counts in the 3-10 keV energy 
range divided by those between 0.7-3 keV, except for 4U 1705-44, 4U 1728-34 and GX 340+0, for which the soft 
band is 1.5-3 keV, 1.8-3 keV and 2.2-3 keV, respectively. 
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Fig. 3. Continued from Fig. [21 



Table 2. Description of the models (M) used to fit the 
spectra in XSPEC notation. 



M Description 

la tbabs*(di skbb+bbodyr ad+gaui /edge i +gau2 ) +gau3 +gau4 
lb tbabs*(diskbb+po+gaui /edgei+gau2)+gau3+gau4 
Ic tbabs*(bbodyrad+po+gaui /edgei+gau2)+gau3+gau4 
Id tbabs* (di skbb+bbodyr ad+po+gaui /edgei+gau2 ) 

+gau3+gau4 
2a tbabs* (diskbb+bbodyrad+gaui /edgei+laor)+gau3+gau4 
2b tbabs* (diskbb+po+gaui /edgei+laor)+gau3+gau4 
2c tbabs* (bbodyrad+po+gaui /edgei+laor)+gau3+gau4 
2d tbabs* (di skbb+bbodyr ad+po+gaui /edgei+laor) 

+gau3+gau4 



power-law component. For example, for the bright source 
GX 349+2, we obtained an index of -0.2 ± 0.4 and an en- 
ergy for the cutoff of I.OIq'j keV. Similarly, for the dim 
:*ource 4U 1728—34 we obtained an energy for the cutoff of 
0.51 ±0.02 keV and the index was unconstrained. Fixing 
the i ndex of the cutoff pow er-law to a more realistic value 
of 1 ( Hasinger et al.lll990[ ) yielded fits with significantly 



and the residuals of the fits are shown in Fig. IC.ll (see 
Appendix [C]) . 

Substituting the blackbody or disc blackbody compo- 
nents in Model la by a power-law component resulted 
in fits with similar or worse xt ^^^ ^^^ the observations 
shown in Table |3] (see Table ^ . Substituting the black- 
body component by a cutoff power-law in Model la re- 
sulted in fits with similar xt- This is expected, since the 
blackbody in Model la accounts for the emission of the 
boundary layer, and could be therefore equally well fit- 
ted by a cutoff power-law or a comptt component, rep- 
resenting saturated comptonization. Finally, substituting 
the disc blackbody component by a cutoff power-law in 
Model la resulted in fits with similar x^, but with unreal- 
istic values for the index and the cutoff energy of the cutoff 



worse xt- 

We performed a comparison of the xt ^'^^ the different 
continua before adding the emission line at ^6.5 keV to 
avoid the EW of the line being affected by the possible 
deficiencies of the fit to the continuum. Whenever two 
different continua yielded a similar xt before inclusion of 
the Fe line, e.g. using a disc blackbody in combination 
"Tvith a blackbody or a cutoff power-law components, the 
breadth of the line was similar in both fits. 

For one of the remaining four observations, corre- 
sponding to XTE J1807— 294, we obtained already an ac- 
ceptable xt of 1-20 for 224 d.o.f. when fitting the contin- 
uum with Model la. However, a better fit was obtained 
substituting the blackbody or disc blackbody components 
by a power law (Model lb or Model Ic, see Table [2]), 
with a final xt of 1-00 (Model lb) and 1.02 (Model Ic) 
for 224 d.o.f. None of the fits showed significant resid- 
ual emission at ^1 or ~6.5 keV. We calculated an up- 
per limit of 17 eV for the EW of a potential Fe line 
with an energy of 6.6 keV and a width of 0.3 keV in 
this observation. For the other AMXP in our sample, 
SAX J1808.4-3658, fitting the spectrum with Models la 
or lb gave significantly worse results {xt of 20.7 and 2.2 
for 224 d.o.f., respectively) than with Model Ic {xt of 1.1 
for 224 d.o.f.). Since XTE J1807-294 was equally well fit- 
ted with Models lb and Ic and SAX J1808.4-3658 was 
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best-fit with Model Ic, we chose the latter to model both 
sources and thus allow a better comparison of the pa- 
rameters. For SAX J1808.4— 3658, we show two fits, cor- 
responding to 2 different spectra: one extracted using the 
full PSF and one after removal of the 2 central columns 
of the PSF (marked with ** in Tables |111S|) . The reason 
is that it was difficult to decide based on the epatplot 
if the original spectrum had some residual pile-up. The 
second spectrum (**) was completely free of pile-up. We 
show the parameters of the best-fit for XTEJ1807-294 
and SAX J1808.4-3658 in Table SI 

Similarly, for 4U 0614+09 we obtained a better xl of 
1.19 (220) with Model Ic, compared to 2.25 (220) with 
Model la and 1.21 (220) with Model lb. We show the 
parameters of the best-fit for this source in Table |4l 

Finally, for Obs 0303250201 of 4U 1636-536 we needed 
three continuum components, blackbody, disc blackbody 
and power-law, in order to obtain an acceptable xt- 
We note that fits of this observation with a simulta- 
neous RXTE observation require as well three compo- 
nents and show a p ower law extending well above 10 keV 
(jPandel et al.ll2008l ). We show the parameters of the best- 
fit for this observation in Table S] 

Next, we examined the asymmetry of the line and eval- 
uated to which extent was the iron emission better ap- 
proximated by a model including relativistic effects. For 
this, we substituted the second Gaussian component at 
~6.5 keV (gau2) in Models la, lb and Ic by a laor com- 
ponent and re-fitted all the spectra with the new models 
(Models 2a, 2b and 2c, see Table [2]). 

We show the parameters of the best-fit with these mod- 
els in Tables [3]|4] and the residuals of the fits in Fig. IC.ll 
The fits with Models 2a, 2b and 2c were again acceptable, 
with xt between 0.8 and 1.2 for a number of d.o.f. between 
161 and 221, i.e. similar to the xt obtained when a simple 
Gaussian model was used to fit the lines. 

We found that fits to continua different to the one 
selected that yielded significantly higher values of xt had 
a clear influence on the line breadth for all the sources 
analysed in this work. 

Narrow absorption features were not visible in any of 
the spectra, as expected for sources which are not at high 
inclination. Similarly, adding an absorption edge to the 
best-fit models shown in Table [3] did not improve sig- 
nificantly the goodness of the fits for any source with 
the exception of GX 349-1-2. For this source, the xt ™" 
proved from 1.15 (216) to 1.07 (214) after adding an edge 
at ^9.2 ±0.1 keV, that we attribute to absorption from 
Fexxvi. The parameters of the Fe line did not change af- 
ter including this additional edge. We note that this edge 
was pre viously detected in a BeppoS AX observation of the 
source (jiaria et al.ll2004l ). 

For sources with significant emission at ^^6.5 keV, the 
parameters of the continuum did not change significantly 
when fitting the former with a Gaussian or a laor com- 
ponent, except for 4U 1728— 34. Previous analyses of this 
source based on C handra and RXTE or BeppoSAX data 
(JD'Ai et al.l 120061 ) showed equally good fits when mod- 



elling the emission at '-~^6.5 keV with a broad line or with 
two absorption edges. Therefore, we regard the XMM- 
Newton fits from 4U 1728— 34 in this work with caution 
and discuss them in the frame of the general properties of 
the iron lines in the following Sections. 

Fig. U shows the ratio between the data and the best- 
fit model at the Fe line energy band for all the sources 
with significant emission when the Fe line is not included. 
The fit of the line with a Gaussian profile is shown in red. 



Table 3. Best-fits to the EPIC pn spectra for all the observations (except for XTE J1807-294, SAX J1808.4-3658, 4U 0614+09 and 
4U1636-536 Obs 0303250201, which are shown in Table S]) using Model la (tbabs*(bbrad+dbb+gaui/edgei+gau2)+gau3+gau4) or Model 2a 
(tbabs*(bbrad+dbb+gaui/edgei+laor)+gau3+gau4). gaui/edgci models the emission feature at ~1 keV with a Gaussian component or an edge respectively. 
fcTbb and fcTdbb are the temperatures of the blackbody and disc blackbody components respectively, febb, fcdbb and fcgau are the normalizations of the blackbody, 
disc blackbody and gaui components in units of (Rin/Dio)^, (R-m/Dio)^ cos9 and 10""' ph cm~^ s~^ respectively. E'gau and a represent the energy and width 
of the Gaussian feature. Eedge is the edge energy and r its optical depth. Fg is the unabsorbed 2-10 keV total flux in units of 10~^erg cm~^ s"'^. We show two 
rows per observation corresponding to fits with Models la and 2a, except for observations for which no significant emission was found at the 6 keV range, since in this case 
Models la and 2a are the same model, namely tbabs*(bbrad+dbb+gaui/edgei)+gau3+gau4. We do not show the values of gau3 and gau4, since they are fixed calibration 
features. The parametres of gau2 and laor components are shown in Tablc[5] The spectrum energy range is 1.5-10, 1.8-10 and 2.2-10 keV for 4U 1705—44, 4U 1728—34 and 
GX 340+0, respectively, and 0.7-10 keV otherwise (see text). 



Target 



Observation 



tbabs 



dbb 



fcTdbb 
(keV) 



bbrad 

feTbb 

(keV) 



fcbb 



gaui 



Fa xi;(d.o.f.) Model 



Jgau 

(keV) 



a 

(keV) 



EW 
(eV) 



4U 1636-536 


0500350301 


2.38 + 0.04 


0.80 + 0.01 


177 + 9 


1.86 + 0.02 


8.24 + 0.3 


1.0 (f) 


0.31 + 0.02 


252 + 30 


63 + 7.5 


1.34 


0.71 


218) 


la 






2.38 + 0.03 


0.80 + 0.01 


177 + 8 


1.86 + 0.02 


8.2 + 0.3 


1.0 (f) 


0.31 + 0.02 


2531?^ 


62 tl 


1.34 


0.70 


217) 


2a 




0500350401 


2.49 + 0.04 


0.84 + 0.01 


185 j:4 


1.80 + 0.02 


12.3 + 0.5 


1.06 + 0.06 


0.23 1°;J^ 


159 ttf 


35 til 


1.82 


0.72 


217) 


la 






9 en +0.09 


0.83+0.01 


196 + 7 


1.80 + 0.01 


12.7 + 0.4 


1-08 tin 


0.20 IH^ 


123 IJf 


27 tf 


1.82 


0.80 


215) 


2a 


GX 340+0 


0505950101 


65.0 + 1.9 


0.96 IH^ 


yrry +235 

'01 -206 


1.64 + 0.03 


14711^ 


- 


- 


- 


- 


15.3 


0.86 


163) 


la 






65.2 IJ;? 


0.95 + 0.05 


787 till 


1.64 + 0.02 


150 til 


- 


- 


- 


- 


15.3 


0.88 


161) 


2a 


GX 349+2 


0506110101 


7.0 + 0.1 


0.92 + 0.03 


465 ir. 


1.57 + 0.02 


169 + 10 


1.06 + 0.02 


0.09 + 0.02 


517ilf,l 


31 tl 


12.1 


1.15 


216) 


la 






7.0 + 0.1 


0.93 + 0.04 


443 + 56 


1.58 + 0.02 


164 1?2 


1.06 + 0.02 


0.09 + 0.02 


519 tilt 


31 + 7 


12.1 


1.14 


214) 


2a 


4U 1705-44 


0402300201 


15.8+°:^ 


1.31 + 0.06 


3.6 + 0.6 


2.76 + 0.2 


0.37 + 0.08 


_ 


_ 


_ 


_ 


0.25 


0.99 


186) 


la 






15.8 + 0.3 


1-30 tHI 


q 7+0.6 
•J.' -0.3 


2.74 + 0.2 


0.38 IHt 


- 


- 


- 


- 


0.25 


1.01 


184) 


2a 




0551270201 


15.0 + 0.5 


0.95 IHI 


231 + 52 


1.61 + 0.03 


79 1^ 


- 


- 


- 


- 


6.5 


1.07 


187) 


la 






14 7+0-5 
J^^-' -0.3 


1.01 + 0.06 


184 t|J 


1.64 + 0.03 


71 + 6 


- 


- 


- 


- 


6.5 


1.03 


185) 


2a 


GX9+9 


0090340101 


1.50 + 0.07 


1.02 + 0.06 


US til 


1.76 + 0.06 


36 + 6 


1.0 (f) 


0.20 + 0.07 


213 + 100 


29+14 


4.55 


0.99 


220) 


la/2a 




0090340601 


1.48 + 0.05 


0.96 + 0.03 


187+17 


1.71 + 0.02 


49 + 3 


1 nt^+003 

'-■"0 _o.05 


0.08 IH^ 


76 til 


lltl° 


5.03 


0.89 


219) 


la/2a 


4U 1728-34 


0149810101 


27.0 + 0.5 


1.9 + 0.3 


3.511- 


3.8 +?1 


0.3 ir2 


- 


- 


- 


- 


0.97 


1.02 


177) 


la 






27.6 + 0.4 


1 48+''-" 
^■^° -0.10 


7.4 IS;^ 


9 "^fi+Oll 
^.''"-0.06 


1.9 + 0.2 


- 


- 


- 


- 


0.97 


1.07 


175) 


2a 


4U 1735-44 


0090340201 


1 sc^+o-io 


0.93 + 0.02 


intll 


1.89 + 0.02 


33 + 2 


1 nA+002 


n ir +0.06 


218 +if 


36 +f 


4.60 


0.68 


216) 


la 






-1 or +0.08 
i.»0_o.o5 


0.93 + 0.02 


172 til 


1.89 + 0.02 


33 + 1 


IMtVol 


0.15 + 0.04 


220 til 


SQtll 


4.60 


0.78 


214) 


2a 


SerX-1 


0084020401 


q f.7+0.08 
•3.01^ _0.06 


0.86 + 0.02 


292 + 24 


1.60 + 0.02 


53 + 3 


1-06 tlf. 


p. iq+0.04 

U.J^'J -0.02 


216 lif 


27 tf 


4.41 


0.98 


216) 


la 






q (!(! +0.08 


0.87 + 0.01 


284 til 


1.61 + 0.01 


52 + 2 


1.06 + 0.03 


T^+O-Oi 
^J-^'J -0.02 


221 ir 


28+ f 


4.41 


0.97 


214) 


2a 




0084020501 


q fin+Oll 


0.91 + 0.02 


233 + 18 


1.64 + 0.03 


33.7 + 2.6 


1 r|y+0.03 

r.ui^ -0.07 


14+0.07 
U.i4 -0.03 


201 tW 


271- 


3.58 


0.93 


216) 


la 






3.591^5 


0.92 + 0.02 


2301^4 


1.65 + 0.03 


33.0 tli 


1 r|y+0.02 


n 14+006 
U.i4 _o 02 


201 +4^ 


271- 


3.58 


0.92 


214) 


2a 




0084020601 


3.68l°:J° 


0.87 + 0.02 


277 til 


1.61 + 0.02 


44.2 + 3.0 


1 06+''°3 
r.UD _o.o6 


n 14+005 
U.i4 _o 03 


203 tlT 


26 +f 


3.97 


1.04 


216) 


la 






3.69l»;°« 


0.87+0.02 


277 til 


1.61 + 0.02 


u.otli 


1 06+0°3 
r.UD _o.05 


n 14+005 
U.i4 _o 03 


201 llf 


251^^ 


3.97 


1.04 


214) 


2a 



AqlX-1 



0303220201 



2.37 



+0.23 
-0.09 



1.10 + 0.03 



25 + 3 



2.25 



+0.07 



3.4 + 0.4 



-, Q-, +0.03 
^.'-'^ -0.09 



0.15 



+0.06 
-0.03 



'" -20 



51: 



1.01 0.97 (221) la/2a 



CD 

CO 

d 

CD 



X 

CO 



edgei 

(keV) 



4U 1543-62 0061140201 



0.71 + 0.04 
0.70 + 0.04 



0.640 + 0.005 
642 +''°°^ 

U.D^Z _o.oo8 



140 + 5 
137 + 5 



1.539 + 0.006 
1.541 + 0.006 



6.991°;!^ 



0.842 



+0.008 



0.20 + 0.02 



6.94 + 0.12 0.841 Iggg^ 0.20 + 0.02 



0.47 0.97 (217) la 
0.47 0.96 (215) 2a 



Table 4. Best-fits to the 0.7-10 keV EPIC pn spectra for XTE J1807-294, SAX J1808.4-3658 and 4U 0614+09 using Model Ic 
(tbabs*(bbrad+po+gaui+gau2)+gau3+gau4) or Model 2c (tbabs*(bbrad+po+gaui+laor)+gau3+gau4) and to the 0.7-10 keV EPIC pn spectra for 4U 1636— 536 
Obs 0303250201 with Model Id (tbabs*(bbrad+dbb+po+gaui+gau2)+gau3+gau4) or Model 2d (tbabs*(bbrad+dbb+po+gaui+laor)+gau3+gau4). The parame- 
ters of the models are the same as in Table |3l except P and fcpi, index and normalization of the power law component, the latter is in units of [ph keV~^ cm~^ 
s"^]. We show one row for each observation whenever no line was found at the 6 keV range, since in this case the two models Ic and 2c or Id and 2d are the 
same model. For SAX J1808.4— 3658 we show results for two spectra, which only differ in the pile-up treatment (see text). 

Target Observation tbabs bbrad po gaui Fg Xi-Cd-O.!.) M 

A^ff" fcTbb fcbb r fcpl Sgau O- fcgau EW 

(10^^ cm~^) (keV) (keV) (keV) (eV) 

XTEJ1807-294 0157960101 5.27 ±0.14 0.69 ±0.02 12 ±2 1.79 ±0.03 0.043 ±0.002 - - - 0.17 1.02 (224) lc/2c 

0.83 0.75 (220) Ic 

0.83 0.71 (218) 2c ^ 

0.80 0.86(220) Ici^ 

0.80 0.84 (218) 2c ^ 



SAXJ1808.4-3658 


0560180601 


1.46 ±0.04 


2.21 ±0.03 


1.12 ±0.08 


2.55 ±0.01 


0.550 ±0.006 






1.43 ±0.04 


2.26 ±0.05 


1 01 +°-^^ 


2.54 ±0.02 


0.546 ±0.006 


SAXJ1808.4-3658 


0560180601** 


1.52 ±0.05 


2.08 ±0.03 


1.4±0.1 


2.59 ±0.02 


0.540 ±0.007 






1.53 ±0.05 


2.08 ±0.03 


1.43 10?5 


2.59 ±0.02 


0.542 ±0.007 



edgei 



gaui 

(keV) (keV) (eV) 



0303250201 


2.77 ±0.03 


1.78 l°-« 


1.0 «i 




^•"J^ -0.05 


1.741^8 
dbb 


1.29 l°;f. 






kTdhb 


fcdbb 






(keV) 








0.42 ±0.04 


63 tg 






0.58l»;^2 


29.9^^5 



td 



(keV) S. 



4U0614±09 0111040101 1.89±0.07 0.64l!5;;|!^ 171^° 2.16±0.02 0.275±0.008 0.856±0.007 0.20±0.02 - - 0.58 0.82(217) Ic S^ 

1.87 ±0.07 0.61 ±0.04 21 +i° 2.15 ±0.02 0.272 ±0.008 0.856 ±0.007 0.20 ±0.02 - - 0.58 0.84 215) 2c b^ 

■ :^ ffi- 



Sgau O- fcgau EW g 



4U 1636-536 0303250201 2.77 ±0.03 1.78+o.t| 1-0 ^^g 1.86 ±0.02 0.2041°;°°^ 1.02 ±0.03 0.221°;°^ 83 ±3 38 ±1 0.76 0.76(216) Id § 



1.8201^:^^1 0.187'I:S"oo4 0.881^;^^ 0.314+^;^^= 1831^^ 691^'^ 0.76 0.78(214) 2d 



P: 



Id 
2d 



X 

CO 



Table 5. Parameters of the Gaussian (Models la, lb, Ic and Id) and laor (Models 2a, 2b, 2c and 2d) components used to model the Fe emission at ^-^6 keV. ii^gau 
and (7 represent the energy and width of the Gaussian feature. The parameters of the laor component E(aor, /?, ri„, i and kiaor account for energy, emissivity 
index, inner radius (in units of r^), inclination and normalization. j3 is constrained to be < 5. The inner radius is constrained to be > 1.235 Tg (GMj (?\ The 
inclination is constrained to be < 70°, since none of the sources is a known dipper or eclipsing system. The outer radius has been fixed to 400 Vg for all sources, 
since it is not well constrained. The normalization is in units of (10^'* ph cm^^ s^^) for both components. Parameters that are completely unconstrained when 
the errors are calculated have been fixed at their best value and are indicated in the table with "(f)" . A star (*) means that the error has reached the maximum 
or minimum allowed value of the parameter. For SAX J1808.4— 3658 we show results for two spectra, which only differ in the pile-up treatment (see text). 



Target 


Observation 


gau2 

-C/gau 

(kcV) 


a 

(keV) 


n^gau 


EW 
(eV) 


laor 

Eiaor 

(keV) 


/? 


r™ 


i 




^laor 


EW 

(eV) 




4U 0614+09 


0111040101 


6.79 + 0.16 


0.94 l°i 


81^ 


185l^J« 


6.50 + 0.10 


2.1 + 0.3 


4n+5.2 

^•U_3.9* 


7Q+0* 


7+^ 


160 tS 


4U 1543-62 


0061140201 


fi 77+0.14 
0-' ' -0.16 


0.33 + 0.13 


0.971°;^^ 


23tji 


0-OD -0.10 


3.6 l^:f 


64l^5l 


631^; 


1.21 + 0.35 


32 + 9 


4U 1636-536 


0303250201 
0500350301 
0500350401 


6.78 + 0.01 
6.97l'^*i6 

r. qy+O.OO. 

O-y' -0.08 


1.15 + 0.17 
0.6 + 0.3 

u.ai _o.26 


14.6 + 9 
3 + 2 

8.3l^:« 


210 + 129 
28 + 19 
5911^ 


6.40 1".^^ 
6.97 (f) 

6.95t»:?r 


2.3 + 0.2 

2.2t\t 
4.8 (f) 


<6.25 

<64 

205.8 (f) 


70+0* 

50 If* 
14(f) 


9 + 1 

q t;+3.3 
0.9 + 0.5 


130+14 
36+?^ 
6 + 3 


GX 340+0 


0505950101 


6.72 + 0.06 


17+0-" 

^•^^ -0.08 


22.0 l^ij* 


17 1^ 


6.82 +H^ 


2.56 (f) 


>16 


<34 


19.7 1^:« 


15 + 5 


GX 349+2 


0506110101 


6.72 + 0.06 


0.32 1°;J^ 


96 i?^ 


82 1- 


5 04 +0.03* 
"■^^ -0.22 


1 +0.7 
^■^ -1.0 


11 -I 


17 + 9 


104121 


96 +S 


4U 1705-44 


0402300201 
0551270201 


6.53 + 0.07 
6.56 + 0.05 


0.34 tH? 
0.44 + 0.06 


1.4 + 0.4 
61 +J'' 


57+16 
92 tf 


6.51 + 0.08 

a 4c +0.05 
"•40 -0.05* 


<3.1 

1 « +0-1 
1^8-0.2 


<110 
<9.1 


45 + 12 

65 +f 


1.4 + 0.4 
83l«3 


58 + 16 
13513,? 


SerX-1 


0084020401 
0084020501 
0084020601 


6.61 + 0.05 
6.58 + 0.07 

O-Oy -0.09 


n 22 +"•"** 

0.25+0.08 

0.27+°;!;; 


22 1« 
15.7 1« 

18-7+5:6 


52+it 
50+1^ 
51+?5 


6.79+0:1 

6.6 + 0.1 

f5 ^0+0.24* 

O-'J -0.23 


HA +0-9 

2^84_o.8 
2.03 (f) 

2.6t?:f 


13.3 i^^ 
11.15 (f) 

11 +26 


<33 

28 + 9 
<25 


25.2 + 4.8 
17.5 + 4 

19.4 i«:^ 


63 + 12 

59 + 13.5 

57+18 


4U 1728-34 


0149810101 


6.72 + 0.10 


0.98+;j;« 


21+.^' 


189 l^f 


a CO +0.11 
"•O-i -0.06 


2.6 + 0.4 


20 1| 


70+0* 


7 7+1.6 
'•' -1.2 


80 122 


4U 1735-44 


0090340201 


6.74 + 0.10 


0.38 + 0.11 


22 + 7 


46 + 14 


6.73 in^ 


<2.94 


5.41 (f) 


44 ir. 


20 + 6 


43 + 13 


SAXJ1808.4-3658 


0560180601 
0560180601** 


6.94 inr 

6.77 + 0.13 


0.26l'^;J^ 
0.39lg;lt 


1 c+0.6 
-■-•0-0.4 

1.9 + 0.7 


23+^ 
30 + 11 


6^50 1°:^^, 
6.96i^:°r 


q +2* 

•J -0.8 

c, n+0* 

^•U-2.6 


6 •^+26 

00 +60 

■J ■J -18 


46 IL 
26 tf 


1 +1.2 
1.8 + 0.6 


CI +20 
OJ^-16 

31 + 10 
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aq 

CD 



CO 
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CD 



CO 
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4. Properties of the sample 

4.1. Fe line emission 

Examination of the spectral fits (Tables [3] to [5] and 
Figs. IC.lHC.4l) shows that we are able to successfully 
model the emission in the Fe K band with a simple 
Gaussian component. Fits of the excess at the Fe K band 
with the more complicated laor component do not im- 
prove significantly the goodness of the fit for any of the 
observations presented here (see Table IH) . 

Considering the whole sample, the Fe line has a 
weighted average energy of 6.67 ±0.02 keV, a width of 
0.33 ± 0.02 keV and an EW of 42 ± 3 eV. The statistical 
distribution of these values is shown in Fig. \5\ The width 
and the EW have a well defined "Gaussian-like" distri- 
bution around the weighted average value. The outliers of 
the histograms in Fig. [S] have values with large errors and 
therefore do not contribute significantly to the weighted 
average. The energy distribution peaks at ~6.7 keV, in 
agreement with the value of the weighted average. This 
is consistent with emission from Fexxv. Clearly, the dis- 
tribution has values consistent with emission from highly 
ionised species of iron, from Fexxii to Fexxvi, and it 
never shows a value consistent with neutral iron. 



MU 
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200 


: 


1 


K. 
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100 


■ l\ 


1 ^ \ 


'f 
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;i 


- 









i 
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1 : 



10*^ 10^ 

Total number of counts 

Fig. 6. EW of the Fe line when fitted with a Gaussian 
component versus the total number of counts of each spec- 
trum. 



In order to find out the reason for the large errors in 
the lines with an EW above 30 eV, we plotted the EW of 
the lines as a function of the total number of counts in the 
spectrum (see Fig. [H]). The lines which show large errors 
are not associated to the spectra with less statistics, in- 
dicating that the detection of the lines within our sample 
is not limited by statistics. Therefore, we examined in de- 
tail the three cases for which the errors on the parameters 
of the line are very large (corresponding to 4U 0614-1-09, 
4U 1636-536 Obs 0303250201 and 4U 1728-34) to deter- 
mine the reason for such large errors. As indicated in 
Sect. 13.21 fits to 4U 1728—34 spectra of Chandra and 



RXTE or BeppoSAX (ID'Ai et al.ll2006[) showed equally 
good fits when modelling the emission at ~6.5 keV with 
a broad line or with two absorption edges. Including one 
edge in the XMM-Newton spectrum at 7.6 keV and no Fe 
line yields a xl of 1.32 (178). Although the xl is higher 
than when including the Fe line, the residuals at the Fe 
band are reduced considerably. The fits to the XMM- 
Newton spectrum have significantly different continuum 
parameters when fitting the residuals at the Fe band with 
Models la or lb (see Table |3]), in contrast to all other ob- 
servations. This indicates that after inclusion of the line, 
the parameters of the fit change significantly and thus the 
line fit may not be realistic. Interestingly, 4U 0614+09 and 
4U 1728—34 show the lowest luminosities of the full sam- 
ple. Therefore, it is plausible that such observations have a 
significantly different, more complex, spectrum compared 
to higher luminosity observations. For comparison, the 
MOS spectra of the observation of 4U 0614-)-09 analysed in 
this work were previously fitted with an absorbe d power- 



law a nd an emission feature at 0.65 keV (Mendez et al 
20021 ) and did not require a line at ~6.5 keV. In the case 



of 4U 1636-536 Obs 0303250201, although the luminosity 
is higher than for other sources in our sample fitted with 
the "standard" Model la/ lb, we needed three continuum 
components to obtain an acceptable fit. This points again 
to a complex spectrum. Fig. [7] shows the residuals of the 
fit to the continuum before including the Fe line for the 
three observations discussed. The residuals at the Fe en- 
ergy band are very different compared to those shown in 
Fig. m indicating that the lines fitted for these observa- 
tions may not be realistic and that such fits may need e.g. 
the inclusion of absorption edges. For these reasons we 
regard the properties of these lines with caution. 

It is interesting to compare the properties of the Fe Ka 
line with general properties of the sources such as luminos- 
ity or temperatures of the blackbody and disc blackbody 
components, with the aim to constrain the origin of the 
line emission region. For example, we may expect a corre- 
lation between the source accretion rate and the breadth 
of the line, if the latter is caused by relativistic effects. 
Fig. [5] shows the properties of the Fe Ka line (width and 
EW) as a function of the luminosity calculated in Table [71 

We do not find any clear correlation between the en- 
ergy centroid, the width or the EW of the Fe line with 
the source luminosity (see Fig. [5]) . The source luminosi- 
ties range over 2 orders of magnitude, covering the full 
range of Eddington ratios from 0.003 to 1.26, while the 
line centroid, width and EW do not show any systematic 
trend. Similary, we do not find any clear correlation be- 
tween the EW of the Fe line and the temperature or flux 
of the blackbody or disc blackbody components. 



4.2. Continuum emission 

We examined the continuum properties of the sample for 
all the sources shown in Table |3l 
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Table 6. xt of spectral fits to (a) diskbb+bbodyrad, (b) diskbb+po and (c) bbodyrad+po components before inclusion 
of the iron line component. Columns (d) and (e) show the xt of the fit after including a Gaussian (d) and a laor (e) 
component to the best-fit model. 



Source 



Observation 
ID 



c 



4U 0614+09 


0111040101 


2.25 


220) 


1.21 


220) 


1.19 


220) 


0.83 (217) 


0.85 (216) 


4U 1543-62 


0061140201 


1.11 


220) 


12.5 


224) 


4.28 


220) 


0.97 (217) 


0.96 (215) 


4U 1636-536 


0303250201 


6.46 


224) 


32.9 


224) 


48.7 


224) 


0.76 (216) 


0.78 (214) 




0500350301 


0.76 


220) 


4.14 


220) 


4.15 


220) 


0.71 (217) 


0.70 (215) 




0500350401 


0.82 


220) 


4.84 


224) 


6.00 


224) 


0.72 (217) 


0.80 (215) 


GX 340+0 


0505950101 


1.15 


166) 


2.04 


166) 


1.83 


166) 


0.86 (163) 


0.88 (161) 


GX 349+2 


0506110101 


1.97 


219) 


3.10 


219) 


3.75 


219) 


1.15 (216) 


1.14 (214) 


4U 1705-44 


0402300201 


1.34 


189) 


1.56 


189) 


1.58 


189) 


1.01 (186) 


1.03 (184) 




0551270201 


3.05 


189) 


5.88 


189) 


4.10 


189) 


1.07 (187) 


1.03 (185) 


GX9+9 


0090340101 


0.99 


219) 


3.58 


221) 


1.26 


219) 


- 


- 




0090340601 


0.91 


219) 


1.27 


221) 


2.13 


221) 


- 


- 


4U 1728-34 


0149810101 


1.55 


180) 


1.82 


180) 


1.80 


180) 


1.02 (177) 


1.07 (175) 


4U 1735-44 


0090340201 


0.99 


219) 


1.82 


219) 


1.98 


219) 


0.68 (216) 


0.78 (214) 


SerX-1 


0084020401 


1.47 


219) 


2.11 


219) 


3.04 


219) 


0.98 (216) 


0.97 (214) 




0084020501 


1.24 


219) 


2.07 


219) 


3.04 


219) 


0.93 (216) 


0.92 (214) 




0084020601 


1.31 


219) 


1.95 


219) 


2.92 


219) 


1.04 (216) 


1.04 (214) 


AqlX-1 


0303220201 


0.97 


221) 


1.22 


221) 


1.35 


221) 


- 


- 


XTEJ1807-294 


0157960101 


1.20 


224) 


1.00 


224) 


1.02 


224) 


- 


- 


SAXJ1808.4-3658 


0560180601 


30.2 


224) 


1.30 


223) 


0.98 


223) 


0.75 (220) 


0.71 (218) 




056018060** 


33.4 


224) 


1.31 


223) 


1.00 


223) 


0.86 (220) 


0.84 (218) 



Figure IHl shows the variation of the blackbody lumi- 
nosity in the band 0.5-30 keV as a function of the total 
luminosity in the same band. In this figure, we excluded 
the two AMSP XTE J1807-294 and SAX J1808.4-3658, 
4U 0614+09 and Obs 0303250201 of 4U 1636-536. The 
reason was that these observations were fitted with a 
continuum different than the standard one of black- 
body and disc blackbody components and are therefore 
indicative of different properties. We further excluded 
the two dim observations corresponding to 4U 1705— 44 
(Obs 0402300201) and 4U 1728-34, since for such dim 
sources Compton scattering is expected to play an im- 
portant role and was not included in our model. 

The blackbody luminosity is very near the value of 
50% of total luminosity expected from simple energy con- 
siderations for all the sources of the sample. Interestingly, 
the temperature of the blackbody component, fcTbb, is 
nearly independent upon the global luminosity at lumi- 
nosities ^2xl0'^^erg s~^ (Fig. IH]). Ignoring the point with 
the highest temperature since it has associated a signif- 
icantly larger error than the other points, fcTbb falls by 
~45% while the total luminosity increases by more than 
two orders of magnitude and the blackbody radius by more 
than a factor of 40. 



5. Discussion 

We performed a systematic analysis of 26 XMM-Newton 
observations corresponding to all the NS LMXBs observed 



with EPIC pn Timing mode and publicly available up until 
the 30th of September 2009 to establish the characteristics 
of the Fc K line emission in these objects. 

In seven observations we did not detect the source sig- 
nificantly. For the remaining nineteen observations we ex- 
tracted one spectrum per observation. We payed special 
attention to the effects of pile-up and background subtrac- 
tion. 



5.1. Fe line emission 

We detected Fe Ka line emission in 80% of the observa- 
tions where the source was significantly detected. 

The energy of the line has values between 6.53 and 
6.97 keV, consistent with highly ionised iron, from Fexxii 
to Fexxvi. The width has values between 0.17 and 
1.15 keV and the EW between 17 and 190 eV. Only in 
four cases out of the fifteen lines reported in this work did 
we find a width near or above 1 keV and an EW above 
100 eV. Interestingly, there is a gap in the distributions of 
a and EW between these high values and the rest of the 
sample. Lines with high values of a or EW ha ve also asso- 
ciated large errors. IVauehan fc UttlevI (|2008l ) interpreted 
a similar finding in the context of detection of narrow 
lines from AGN as the lines with large EWs being most 
likely false detections. In the cases studied here the lines 
with widths near 1 keV are also in the limit of detectabil- 
ity of 3(7 and the large errors most likely point to an in- 
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Fig. 4. Ratio of the data to the continuum model for ah the neutron star LMXBs analysed in this work for which 
significant Fe K emission is detected. 
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Table 7. 2-10 keV luminosities in units of 10'^^ erg s ^ and corresponding Eddington ratios (L2-10 kev/^Eddington) 
for all the observations in the sample. 



Source 


Observation 


d 


L36 


Rfidd 




ID 


[kpc] 


erg s"^ 




4U 0614+09 


0111040101 


3.2 + 0.5 fKuulkers et al. 2009") 


0.73 


0.004 


4U 1543-62 


0061140201 


7 (Wane & Chakrabartv 2004) 


2.87 


0.016 


4U 1636-536 


0303250201 


6.0+0.1 fGallowav et al. 2008) 


3.41 


0.019 




0500350301 




6.01 


0.033 




0500350401 




8.17 


0.045 


GX 340+0 


0505950101 


11 fGrimm et al. 2002) 


230.72 


1.264 


GX 349+2 


0506110101 
0402300201 


9.2 fGrimm et al. 2002) 


127.64 
1.05 


0.701 


4U 1705-44 


5.8+0.2 fGallowav et al. 2008) 


0.006 




0551270201 




27.25 


0.150 


GX9+9 


0090340101 
0090340601 


4.4 fGrimm et al. 2002) 


10.98 
12.14 


0.060 






0.067 


4U 1728-34 


0149810101 


4.0 + 0.4 (Galloway et al. 2008) 


0.48 


0.003 


4U 1735-44 


0090340201 


6.6 + 1.0 (Galloway et al. 2008) 


24.97 


0.137 


SerX-1 


0084020401 


7.7 + 0.9 (Galloway et al. 2008) 


32.59 


0.179 




0084020501 




26.45 


0.145 




0084020601 




29.33 


0.161 


AqlX-1 


0303220201 


3.9 + 0.7 (Gallowav et al. 2008) 


1.91 


0.010 


XTEJ1807-294 


0157960101 


8 (assumed) 


1.36 


0.007 


SAXJ1808.4-3658 


0560180601 


2.77 + 0.11 fGallowav et al. 2008) 


0.79 


0.004 
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Fig. 5. Histogram of the parameters of the Fe line when fitted with a Gaussian component for the full sample considered 
in this work. The white rectangles marked with an arrow correspond to the lines which are detected below a significance 
of 3 cr (see Fig. [5]). The width (middle panel) and the EW (right panel) have a well defined "Gaussian-like" distribution 
around the weighted average value, while the outliers have large errors associated to their a and EW values. 
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Fig. 7. Residuals in units of standard deviation from the best-fit model to the data for 4U 0614+09, 4U 1636— 536 
Obs 0303250201 and 4U 1728-34 before including the Fe line. 
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Fig. 8. Energy, EW and a of the Fe line versus luminosity of the source in the 2-10 keV energy band. The luminosity 
has been calculated based on the distances reported in Table [T] 



m 100 r 




10 -00 

Luminosityj^i [lO^'ergs/s] 



10 -00 

Luminosityj^i [lO^'ergs/s] 




Rodius,, [km] 



Fig. 9. Left panel: Variation of blackbody luminosity in the band 0.5-30 keV with total luminosity in the same band. 
The dashed line shows the 50% of the total luminosity. Middle panel: Temperature of the blackbody component versus 
the total luminosity in the 0.5-30 keV energy band. Right panel: Temperature of the blackbody component versus the 
radius of the blackbody emission. The radius has been inferred from the normalisation of the blackbody component. 
See Sect. 14.21 for the observations included in this figure. 



appropriate modelling of the continuum for 4U 0614+09, 
4U 1728-34, and 4U 1636-536. 

Considering the whole sample, the Fe line has a 
weighted average energy of 6.67 ±0.02 keV, a width (cr) 
of 0.33 ± 0.02 keV and an EW of 42 ± 3 eV. 

Recently, broad skewed Fe K emission lines from 
the disc were reported from LMXBs containing a NS 



(e.g . Bhattachar vva fc Strohmaveiil2007t llaria et al.ll2009t 
D^i et al. 2009). All the XMM-Newton observations for 
which a broad skewed Fe line has been reported were in- 
cluded in our sample (see Sect. S]). However, in contrast 
to previous analyses, we do not need to invoke relativis- 
tic effects to explain its width, which could be due to 
mechanisms such as Compton broadening. The lines are 
equally well fitted with the relativistic laor component 
or with a simpler Gaussian component. Further, the pro- 
files shown in Fig. 3] do not show an asymmetric shape, 
as expected if the lines are emitted close to the NS and 
shaped by relativistic effects. The line profile is instead 
symmetric, simila.rly to the one found in dipping sources 
( Diaz Trigo et al.l2006r) . The major difference between the 
analysis presented in this work and previous works is a 
careful treatment of pile- up effects, common in the obser- 
vations of bright LMXBs, and a different modelling of the 
continuum in some cases, which has a strong effect in the 
EW of the lines. 



As shown by iBrandt fc MattI ( 1994 ). in a scenario 
where the iron line is generated at the inner disc, we would 
expect the centroid and the width of the iron line to de- 
crease as the accretion rate diminishes and the disc re- 
cedes. The behaviour of the equivalent width is somewhat 
more complex: the line equivalent width is expected to first 
increase with ^ up to a cer tain value at whi c h it s tarts to 
decrease (see Fig. 5b from iBrandt fc MattI (|l994l )l. This 
is due to the emission properties of iron atoms as a func- 
tion of the ionisation parameter (Matt et al. 1994). When 
the most abundant ions are Fexxiv-Fexxvi, line emis- 
sion is more intense than in the neutral case, due to the 
greater fluorescent yield and the smaller photoabsorption 
cross-section at the line energy. However, if ^ increases 
further, the ion distribution becomes dominated by fully 
stripped atoms and the line emission decreases. In the 
observations presented in this paper, we do not observe 
a correlation between the centroid and the width of the 
iron line with luminosity, in agreement with a similar sys 



tematic analysis based on ASCA observations (jAsai et al. 



20001 ). In contrast, the three cases for which a large width 



and equivalent width are measured occur at relatively low 
luminosities and are actually associated to large errors. 
Th is seems to be at odd s with the behaviour predicted 
by IBrandt fc MattI ( 1994 ) for lines generated at the disc 
and suggests instead that the physical condition of the 
line-emitting region is rather similar among the LMXBs. 
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An interpretation of the line in the context of disc 
emission would imply that the line is produced far from 
the neutron star in order to explain its symmetric shape. 
However, in such a case it is difficult to justify the varia- 
tions of the line as the disc inner radius changes. 

Alternatively, the origin of the line broadening could 
be Compton scattering in a corona. Based on a system- 
atic an alysis of twenty NS LMXBs with ASCA, lAsai et al. 



(120001) concluded that the iron lines were likely produced 
through the radiative recombination of photoionised 
plasma, being the line width a result of the combination 
of line blending, Doppler broad ening and Compton scat - 
tering. Combining such effects, iKallman fc White! (J1989I) 
estimated that EWs of up to 100 eV were attainable from 
an ADC using standard parameters. This in agreement 
with the values obtained in this work. 

Finally, the line could originate in a partially ionised 
wind as a result of illumination by the central source con- 
tinuum photons and broadened by electron downscatter- 
ing in the wind environme nt (e.g., iLaurent &: TitarchukI 
2007tlTitarchuk et al.ll2009l) . The line profiles obtained via 
such mechanism are again asymmetric, similar to the ones 
produced by relativistic effects. Since the line profiles ob- 
tained here are highly symmetric, we do not attribute the 
line origin to the outfiow based on these observations. 

5.2. Continuum emission 

The X-ray continua for 80% of the observations were well 
fitted by a model consisting of a blackbody and disc- 
blackbody components absorbed by neutral material. For 
15% of the observations, a model consisting of a black- 
body and power-law components absorbed by neutral ma- 
terial was preferred to the combination of blackbody and 
disc-blackbody components. For one observation a fit with 
three components, namely, blackbody, disc-blackbody and 
power- law was required in order to get an acceptable xt- 
For the fits with blackbody and disc-blackbody com- 
ponents the temperature of the blackbody component had 
values between 1.5 and 2.8 keV, except for 4U 1728— 34, 
for which a temperature of 3.8 ~ti'i keV was found. The 
temperature of the disc blackbody component had lower 
values between 0.6 and 1.3 keV, except for 4U 1728— 34, 
for which a temperature of 1.9 ±0.3 keV was found. 
This is consistent with the interpretation of the black- 
body component as the boundary layer between the neu- 
tron star and the inner di sc, which shows smaller radii 
and h igher temperatures. IChurch fc Baluciiiska-Church 
(|200ll) . based on an ASCA survey, found a broad de- 
pendence of the blackbody luminosity on the total lu- 
minosity not previously known for LMXBs in gen- 
eral, but while they found that the blackbody lumi- 
nosity was falling more rapidly than the total lumi- 
nosity as the mass accretion rate decreases, we find 
a very tight correlation where the blackbody luminos- 
ity is very near the value of 50% of total luminosity 
at all luminosities. This may be due to the different 



continuum m odel used, blackbody and cutoff po wer-law 
components ([Church fc Balucinska-ChurchI 120011 ) versus 
blackbody and disc blackbody components (th is work). 
Similarly to lChurch fc Balucihska-ChurchI ( 200lh . we find 
that the major change in blackbody emission results pri- 
marily from changes in the emitting area, not the tempera- 
ture. The independence of the temperature of this compo- 
nent on the global value of Eddington ratio lends support 
to the theoretical suggestion that the boundary layer is 
r adiation press ure supported (jlnoeamov fc Sunvaevll 19991 : 
iRevnivtsev fc Gilfanov 200(1 . 

For the sources fitted with blackbody and power-law 
components, the temperature of the blackbody component 
had values between 0.6 and 2.2 keV and the index of the 
power-law was between 1.8 and 2.6. 

We detected an excess of emission at ~1 keV in all 
the observations where we could use the spectrum down 
to 0.7 keV, except for the two accreting millisecond pul- 
sars XTEJ1807-294 and SAX J1808.4-3658. We mod- 
elled such emission with a Gaussian component centred 
between 1.0 and 1.15 keV in 12 observations and as 
an edge between 0.84 and 0.86 keV in 2 observations. 
Provided the feature has an astrophysical origin, its ap- 
pearance always at similar energy points to line emission, 
since a soft component due to a e.g. blackbody emission is 
expected to change its temperature for different sources. 
Its energy is consistent with emission from highly ionised 
Fe, from Fexx to Fexxiv, and/or Nex. We did not find a 
clear correlation between the width of the line emission at 
~1 keV and the Fe Ka band. However, the fact that the 
same feature is not observed in simultaneous EPIC and 



RGS exposures (e.g. iBoirin et al.ll2005l ) indicates that the 
feature is likely blended emission from various ions and 
additionally broadened by the same mechanism that the 
Fe Ka line. 

In this work, we show that effects of pile-up can 
significantly affect the width, and therefore the physi- 
cal interpretation, of the iron line in NS X-ray bina- 
ries. Furthermore, parallel works on XMM-Newton and 
Suzaku observations of the black hole X-ray binary 
GX 339-4 demonstrated that, taking into account pile- 
up effects, the breadt h of the li ne was consistent with a 
truncation of the disc (JDone fc D iaz Trigo 2009), or with 
a black hole with moderate spin (lYamada et al.l 2pp9[). i n 
contrast to previous claims (e.g. iMiller et al.ll2006ll2008l ). 
Therefore, we urge caution in using piled-up data for de- 
tailed spectral analysis. 
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Appendix A: Examples of pile-up removal 



Fig. lA.ll and IA.2I show the effects of pile-up in the 
epatplot in the bright source GX 349+2 and the dim 
source 4U 1728—34, respectively. 

For GX 349+2 we observe a significant energy- 
dependent deviation between the data and the model > 4 
keV when the full PSF region is used (Fig. lA.ll upper 
left). 

As we remove the inner columns from the PSF the ef- 
fects of pile-up are mitigated. After removal of 6 columns, 
there is still a small deviation above 9 keV. After removal 
of 8 columns, there is not an energy-dependent deviation 
anymore (Fig. lA.ll right panel). 

For 4U 1728—34 we do not observe any energy- 
dependent deviation between the data and the model 
when the full PSF region is used (Fig. IA.2[ left panel). 
Similarly, when we excise the core of the PSF the pattern 
ratios do not change with respect to the case where the 
core of the PSF is included. This is an indication that 
pile-up is not present in this source. 

Appendix B: Comparison with previous results 

A number of the sources analysed in this pa- 
per were previously published by o ther author s (e.g . 



Bhattacharvva fc Strohmaven 120071 : llaria et al.l l2009t 



D'Ai et a l. 2009) . In what follows we compare the proper- 
ties of the Fe Ka line found in previous analyses with the 
ones inferred in this paper. 

For most of the sources we obtained significantly dif- 
ferent parameters for the Fe line when fitted with a laor 
component, compared to previous analyses (we discuss the 
discrepancies in detail below). 

Therefore, in order to obtain a quantitative measure- 
ment of the difference between the Fe line obtained by 
other authors and in this work we performed fits to the 
spectra with the best-fit model shown in Sect. 13.21 but 
fixing all the parameters of the laor component to the 
values obtained in previous analyses except the normali- 
sation, which was left free. The xl of these fits is shown in 
Table IbTTI For comparison, we show in the same table the 
xt obtained when the line is fitted both with a Gaussian 
and a laor component and leaving all the parameters free, 
as in Tables ^SSi We note that we compar ed the results 
in this work with the results obtained by ICackett et al. 
( 2009a) only in Table IB. II and not in the subsections that 
follow, since the latter work is in refereeing process at the 
moment of acceptance of this paper, and their parameters 
may still change in their final version. 

We found that the xt oi the fits to the lines with the 
parameters of the literature yielded higher values of xt 
for eight of the observations analysed in this work (see 
Tabic EB. For the three observations of 4U 1636-536, 



the difference in xt i^ ^^ot statistically significant. Indeed 
the lines fitted in this work have a width a near 1 keV, 
similarly to previous analyses. However, we note that two 
of the three lines are below the level of 3 cr detectabil- 
ity significance in this work and the third one is only 
slightly above such limit (see Sect. U). We also note that 
in the case of 4U 1636-536 Obs 0500350401, we obtained 
a lower, though not significant, xt when using the param- 
eters of the line from previous analyses. This is due to the 
fact that we did not allow in our analysis an inclination 
higher than 70°, since the source shows neither dips nor 
eclipses. In contrast, the value of the inclination in the 
laor compone nt in previou s analyses had a large valu e 



of >81° fe.g.. lPandelet al 



2008 : Cackett et al.ll2009bh . 
For the observations of SerX-1, SAX J1808.4-3658 and 
GX 340+0, the increase of the xt when we use the pa- 
rameters of the lines in the literature is not statistically 
significant, since the xt is already well below 1 in the fits 
in this work. However, we note that since the xt i^ already 
below 1 for fits with the simple Gaussian model there is 
no reason to consider a laor component for these fits. 

B.l. SerX-1 



iBhattacharvva fc: Strohmaveij (2007) first reported the 
presence of a skewed iron Ka line with a moder- 
ately extended red wing in a NS LMXB based on the 
three XMM-Newton observations of SerX— 1 that we 
re-analysed in this paper. They fitted the line with a 
laor component and found EWs for the lines between 



86 



-1-9 



and 105 li cV. 



16.19 



-1-8.77 
-3.19 



inner radii between 4.04 tpgg and 
-0.0* ^^^■ 



Tg and an energy centroid of 6.4 



On average we found smaller EWs, between 57+ 18 and 
63 + 12 eV and larger energy centroids, between 6.60 ±0.1 
and 6.79 lg'2 keV, when fitting the line with the same com- 
ponent and EWs between 50l]^2 ^^^ 52li[j^ and an aver- 
age energy centroid of 6.60 ± 0.04 keV when fitting the line 
with a Gaussian component. We also inferred a smaller 
average inclination for the source from the laor fit: be- 
tween <25 and 28 + 9°, compared to between 39.7tJ_^and 



50.2ttl° found bv IBhattacharvva fc Strohmaveij (|2007l ). 



The inner radius was not well constrained when fitting the 
line with a laor component, with a value of 13.3 _g g ig 
for the best constrained observation (Obs 0084020401). 
In our analysis the fit with the laor component was not 
statistically preferred to the one with the Gaussian com- 
ponent. 

We explain the discrepancies in the results as due 
mainly to two differences in the analysis: first, while 



Bhattacharvva fc Strohmaveij (|2007^ considered pile-up to 



be insignificant for these observations, we found that the 
epatplot shows a significant degree of pile-up before we 
remove the central 4 columns of the PSF. Fig. IB. II shows 
the difference in the spectra extracted including the core of 
the P SF (red), as done bv IBhattacharvva fc Strohmavei 
(J2007l ). and excluding the 4 central columns to remove 
pile- up effects (black), as in this paper. The pile- up effects 
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Fig. A.l. Measure of pile-up in the EPIC pn Timing Mode using the SAS task epatplot for the bright source 
GX 349+2. Source events were extracted from a 64" (16 columns) wide box centred on the source position (left 
panel). Middle-left, middle- right and right panels show the same region as in the left panel but after exclusion of the 
neighbouring 4, 6 and 8 columns from the centre of the box. 




Fig. A. 2. Same as Fig. lA. II but for the dim source 4U 1728-34. 



show up as an excess of photons from as low as 3 keV and 
becomes especially prominent above 5 keV in the red spec- 
trum. This, combined with a improper continuum mod- 
elling (see below), may create an excess of photons below 
6 keV which could be misinterpreted as a red wing of the 
Fe line. 

The second difference in the analysis is precisely that 
we used a different continuum model to fit the spectrum. 
As shown in Table [S] we got a significantly worse xt using 
a combination of diskbb and comptt components to fit 
the co ntinuum, as done bv iBhattacharvva fc Strohmaver 
(J2007() . when compared to diskbb and bbody components, 
as used in this paper. 

Therefore, for this particular source, the main contri- 
bu tion to the broad Fe line detec t ed in the first analysis 
by IBhattacharvva fc: Strohmaven (|2007l ) is likely related 
to residuals from the fit to the continuum, showing that 
modelling the spectra with different continua may modify 
significantly the parameters of the Fe line. 



B.2. GX 340+0 



DAi et al.l (|2009[ ) reported the presence of a broad asym- 
metric emission line in the XMM-Newton spectrum of 
GX 340+0. Due to the variability of the source as revealed 
by the light curve (see Fig. [5]), they divided the observation 
in five segments and analysed each segment separately. 
They modelled the spectral continuum with diskbb and 
bbody components and found a significant improvement 
in the fit when using the diskline component to model 
the Fe Ka emission compared to a Gaussian component. 
The line had an average energy of 6.69 + 0.02 keV, in- 
ner radius of 13+3 ig, and EW of 41 + 3 eV. In spite of 
using the same continuum to model the spectrum, we ob- 
tained equally good fits when modelling the Fe line with a 
Gaussian or with a laor component. The line had an av- 
erage energy of 6.82 1013 keV, an inner radius >16 r^, and 
a smaller EW of 15 + 5 eV. We inferred an average incli- 
nation for the source of < 34°, compared to 35+1° found 
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Table B.l. xt of spectral fits for sources for which asymmetric Fe lines showing relativistic effects have been reported 
in the literature. Columns (a) and (b) show the xt of the fit after including a Gaussian and a laor component, 
respectively, to the best-fit model. Column (c) shows the xt of the fit after including a laor component with all the 
parameters fixed to the values obtained by previous authors, except the normalisation, to the best-fit model. 



Source 



Observation 
ID 



xl 



4U 1636-536 



GX 340+0 

GX 349+2 

4U 1705-44 

SerX-1 



SAXJ1808.4-3658 



0303250201 
0500350301 
0500350401 
0505950101 
0506110101 
0402300201 
0551270201 
0084020401 
0084020501 
0084020601 
0560180601 



0.76 (216) 
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Fig. B.l. Ratio of the EPIC pn spectrum of SerX— 1 
(Obs 0084020501) to its best fit continuum model for the 
spectrum free of pile- up used in this paper (black). The 
red points show the ratio of the piled-up spectrum of the 
same observation, obtained when the full PSF is used, to 
the best fit continuum model of the spectrum free of pile- 
up. The pile-up shows up mainly as a significant hardening 
of the spectrum at energies ^5 keV. This adds an 'artifi- 
cial' red-wing to the iron line due to the different curvature 
of the spectrum. 



bv ID'Ai et al.l ( 20091) . When fitted with a Gaussian pro- 
file, the line had a n energy of 6.76 ± 0.02 keV and width 
of 0.24 ±0.02 keV (JD'Ai et al.ll2009l ). consistent with our 



values of 6.72 ± 0.06 keV and width of 0.17 to.lt keV. 

The main difference bewteen our analysis and theirs 
is again the treatment of pile-up. While they considered 
only the initial ^7 ks of the observation (interval 5) to 



be affected by pile-up and excised the 2 central columns 
to remove its effects, we found that the full observation 
was strongly affected by pile-up and removed the central 
8 columns. This is consistent with the high count rate of 
the source, with peaks up to ^1100 counts s~^and a min- 
imum of 650 counts s~^ along the observation, as well as 
the hardness of the spectrum (see Sect. 12. 1( ). Note that the 
light curve in this paper (see Fig. [2]) shows a count rate 



twice as large as the one shown in Fig. 1 from ID'Ai et al 



( 20091 ). A potential expla nation of this diff erence in the 
light curves could be that ID'Ai et al.l ( 20091 ) averaged in- 
tervals of real data with those for which no data were 
recorded, due to telemetry saturation. 

B.3. GX 349+2 



laria et al.l ( 20091) reported the presence of several emis- 
sion features, consistent with the transitions of L-shell 
Fe XXII- Fe xxiii, S xvi, Ar xviii, Ca xix and highly ionised 
Fe in the XMM- Newton spectrum of GX 349+2. While 
the first four features could be fitted equally well with 
Gaussian features or the relativistic diskline component, 
they found that the Fe Ka feature was better fitted using 
the diskline component at 6.76 keV or two diskline 
components at 6.7 and 6.97 keV. The line had an energy 
of 6.80+0.02 keV, width {a) of 0.28tn°il keV and EW of 



-0.04 



49_7 eV when modelled with a Gaussian component. The 
same line had an energy of 6.76±0.02 keV, inner radius 



6.2 



+ 19.1 
-0.2 



Ya and EW of 61±9 eV when modelled with a 



diskline component. 

In spite of using the same continuum to model the 
spectrum, we obtained equally good fits when modelling 
the Fe line with a Gaussian or with a laor component 
and a significantly larger EW in any of the fits, compared 
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to llaria et al.l ( 20091 ) . In this paper, the hne has an energy 
of 6.72±0.06 keV, width (a) of 0.32t|]J^ keV and EW 
of 82^^5 eV when modehed with a Gaussian component 



and energy of 6.94 



+0-03 1^ Y 
0.22 ^^^^ 



inner radius of 11 



5 ^9 



and 



EW of 96^^4 eV when modehed with a laor component. 
We obtained a significantly different i nclination for th e 
source, 17±9°, compared to 41.4^2.1° (|laria et al.ll2009t ). 
The line parameters, apart from the inclination and the 
EW, were not significantly different within the errors. We 
attribute the differences in some prope rties of the l ine to 
the PSF extraction regions used. While llaria et al.l (|2009r) 



excluded only the 2 central pixels of the PSF before spec- 
tral extraction, their Fig. 4 shows an excess of double 
events and a deficit of single events above ~7.5 keV com- 
pared to the predicted model. This is a clear indication of 
residual pile-up in their spectra (see XMM-Newton Users 
Handbook). In contrast we found that rejection of the in- 
ner 8 columns of the PSF was necessary in order to remove 
the pile-up effects completely. 

We found a lower temperature for the diskbb and 
bbody components, 0.92 ±0.03 and 1.5 7±0.02 keV, com - 



pared to 1.05 tHa and 1.7921°:°?^ keV (llaria et al.ll2009h . 



Thus the spectrum extracted in this p aper is sign if icantl y 
softer compared to the analysis by llaria et al.l (|2009t ). 
Again, this is expected if pile-up effects have not been 
completely removed from the spectra in their analysis. 



B.4. 4U 1705-44 

4U 1705-44 was observed twice by XMM-Newton with 
the EPIC pn camera in Timing mode. 

The first observation had a relatively low flux, 
2.5xl0~'^°erg cm~^ s~^, and was at least one order of 
magnitude below the typical count rate where pile-up ef- 
fects become important. In 2008 XMM-Newton observed 
4U 1705— 44 at a significantly higher fiux, 642 counts s"^ 
compared to 28 co unts s~^ in 2006 (see Table [J). 
Di Salvo et al.l ( 20091 ) reported the presence of four emis- 
sion features and one absorption edge in the spectral resid- 
uals after modelling the continuum with a bbody and a 
comptt components, which were identified with emission 
of Sxvi, Ar XVIII, Caxix and Fexxv, and (redshifted) 
absorption from Fexxv. The lines were broad with a be- 
tween 120 and 260 eV. They found a significant improve- 
ment of the fit when the Fe xxv feature was modelled with 
a diskline component compared to a Gaussian compo- 
nent and reported values of 6.66 ±0.01 keV, 14 ±2 r^, 
39 ± 1 and 56 ± 2 eV for the energy, inner radius, incli- 
nation and EW for the fit with a diskline component. 
The absorption edge appeared to be smeared (width of 0.7 
keV) and redshifted with an energy of 8.3 ± 0.1 keV with 
respect to the rest frame energy of 8.83 keV. We found 



values of 6.45 ±0.05 keV, <9 ig, < 34° and 135^2? eV 
for the energy, inner radius, inclination and EW for the 
fit with a laor c omponent, i.e. signifi cantly different to 
those reported bv lOi Salvo et al.l (|2009l ). When fitted with 



a Gaussian component the centroid was 6.56 ± 0.05 keV 
and the EW 92 tf eV. 

The difference between both analyses is related again 
to both the pile-up treatment and the modelling of the 
continuum. We found that it was necessary to excise the 7 
central columns of t he PSF in order to rem ove completely 
pile-up effects, while IPi Salvo et al.l (1 2009) did not remove 
any column. In addition. lPi Salvo et al.. (,2009. ) used a con- 
tinuum of bbody and comptt to fit the spectrum, while we 
used bbody and diskbb. 



B.5. SAX J1808.4- 3658 

SAXJ1808.4-3658 was observed by XMM-Newton in 
2008 and its spectrum modelled with a combina- 
tion of diskbb, bla c kbody and power l a w components 



( Paoitto et aD 120091: IPatruno et all 120091: ICackett et al 



2009al ). The residual emission at the Fe Ka band is mod- 
elled with a d i sklin e profile in all the analyses, although 
Papitto et al.l ( 2009 ^ obtained only a slight improvement 
in the fit when using a diskline profile compared to a 
Gaussian profile. The parameters of the diskline compo- 
nent were in general consistent in the three analyses, with 
some exceptions such as the inclination derived for the 
source, due probably to small differences in the analyses 
such as including or not the RGS in the simultaneous fit or 
considering shghtly d ifferent energy bands . For example. 



Papitto et al 



20091) (jCackett et al.ll2009al) used the 1.4- 
11 keV (1.2-11 keV) band for spectral fitting since they 
found large residuals below 1.4 (1.2) keV, respectively. 

As explained in Sect. 13.21 we extracted for this obser- 
vation two different spectra. Spectrum 1 was extracted 
including the core of the PSF, as done by the authors 
above. Despite using the same extraction region, we found 
a significantly smaller line, with an EW of 51 l^g eV when 
fitting the line with a laor component and 23 _7 eV in the 
fits with a Gaussian component. In contrast, in the previ- 
ous analyses with a laor component, the EW was found 
to be as large as 121 + ?g eV (Papitto et al. 2009). 118 ± 10 



( Cackett et aLll2009al) and 97.7±31.4 eV (jPatruno et al 
20091) . 



This may be related to the different continuum used 
to fit the spectrum. While all the previous analyses of the 
XMM-Newton observation of SAX J1808.4-3658 required 
3 continuum components to achieve an acceptable fit, we 
obtain already a xt of 0-98 (223) with a continuum con- 
sistent of blackbody and power-law components. 

B.6. 4U 1636-536 



Pandel et al.l (|2008l ) found evidence for the presence of 
relativistic lines from iron in different ionization states in 
the three XMM-Newton observations of 4U 1636—536 pre- 
sented here. They reported EWs of 215, 98 and 140 eV for 
the iron Ka line present in Obs 0303250201, 0500350301 
and 0500350401 respectively, which they modelled with a 
diskline component. In contrast we found smaller EWs 
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of 130 ± 14, 36 tit and 6 ± 3 eV when fitting tlie line with 
a laor component and of 210 ± 129, 28 ± 19 and 59 ^{7 eV 
when using a Gaussian component instead. We note that 
when fitting the line with a Gaussian component, the lines 
from Obs 0303250201 and 0500350301 are below the 3 a 
significa nce and should ther efore not considered as de- 



tections. iPandel et al.l (|2008r ) regarded the high inclina- 



tion, > 81°, obtained for Obs 0303250201 and 0500350401 
as unrealistic (given the valu es of the inclination of 36- 
74°determined for this source ICasares et al.l ( 20061) ) and 



interpreted it as an indication that the excess at the Fe Ka 
band was a blend of at least two lines. In contrast, we 
found only a high inclination for Obs 0303250201 and an 
upper value of 70° (as used for all the fits in this sam- 
ple) gave already an acceptable fit with a xt of 0.76 (216) 
for this observation. We obtained equally good fits when 
modelling the Fe line with a Gaussian or with a laor com- 
ponent. 

We attribute the differences in some pro perties of the 



line t o the PSF extraction regions used. iPandel et al 



pe 

i 



( 2008h considered that none of the three XMM- Newton 



observations suffer from pile-up effects based on the count 
rate limit for the pn Timing mode. In contrast we found 
that rejection of the inner 1(3) columns of the PSF was 
necessary in order to remove the pile-up effects completely 
for Obs 0500350301 (0500350401). The residual pile-up ef- 
fects in their spectra are most likely responsible for the 
harder spectrum. They reported a difference of their spec- 
tra with respect to the simultaneous RXTE PGA spectra 
both of flux and slope. The PGA spectra show a flux excess 
of ~30% at 3 keV and -10-15% at 10 keV with respect to 
their EPIG pn spectra. Qualitatively, this is what we ex- 
pect if the XMM-Newton spectra are affected by pile-up: 
photons will be lost so that in average the flux is lower and 
in addition soft photons will be counted as hard photons, 
so that in average the PGA will show a softer spectrum. 

Appendix C: Spectra 

The best-fits of each individual spectrum together with 
the residuals of the fits are shown in Figs. IG.mG.2l 
Figs. IG.3IIG.4I show the unfolded spectra for each obser- 
vation. 
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Fig. C.l. Upper panels: 0.7-10 keV EPIC pn (black) spectra fit with Models la-Id. For each source the best-fit is 
shown (see Tables [3]|4] and text). Lower panels: Residuals in units of standard deviation from the corresponding model. 
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Fig. C.2. Continued from Fig. IC.ll 
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Fig. C.3. 0.7-10 keV EPIC pn (black) spectra fit with their best-fit model (see Tables [3]|3]) shown in flux units. 
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Fig. C.4. Continued from Fig. lC.3l 



